Dissipative dynamics of topological defects in frustrated Heisenberg spin systems 
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We study the dynamics of topological defects of a frustrated spin system displaying spiral order. 
As a starting point we consider the SO(3) nonlinear sigma model to describe long-wavelength 
fluctuations around the noncollinear spiral state. Besides the usual spin-wave magnetic excitations, 
the model allows for topologically non-trivial static solutions of the equations of motion, associated 
with the change of chirality (clockwise or counterclockwise) of the spiral. We consider two types of 
these topological defects, single vortices and vortex-antivortex pairs, and quantize the corresponding 
solutions by generalizing the semiclassical approach to a non-Abelian field theory. The use of the 
collective coordinates allows us to represent the defect as a particle coupled to a bath of harmonic 
oscillators, which can be integrated out employing the Feynman- Vernon path-integral formalism. 
The resulting effective action for the defect indicates that its motion is damped due to the scattering 
by the magnons. We derive a general expression for the damping coefficient of the defect, and 
evaluate its temperature dependence in both cases, for a single vortex and for a vortex-antivortex 
pair. Finally, we consider an application of the model for cuprates, where a spiral state has been 
argued to be realized in the spin-glass regime. By assuming that the defect motion contributes 
to the dissipative dynamics of the charges, we can compare our results with the measured inverse 
mobility in a wide range of temperature. The relatively good agreement between our calculations 
and the experiments confirms the possible relevance of an incommensurate spiral order for lightly 
doped cuprates. 

PACS numbers: 75.10.Nr, 74.25.Fy, 74.72.Dn 



I. INTRODUCTION 

Two-dimensional frustrated Heisenberg spin systems 
with noncollinear or canted order have attracted much 
attention recently. Noncollinear order arises due to frus- 
tration, which may originate from different sources. The 
most common kind of frustration is realized in antifer- 
romagnets on a two-dimensional (or three-dimensional 
stacked) triangular lattice. Prototypes of these geomet- 
rically frustrated magnets are pyrochloresi 2 ^ A sec- 
ond source of frustration may be a competition between 
nearest-neighbor and further-neighbor exchange interac- 
tions between spins. Typical examples are helimagncts, 
where a magnetic spiral is formed along a certain direc- 
tion of the latticed A third kind of frustration may occur 
by chemical doping of a magnetically ordered system. In 
this case, the spin-current of the itinerant doped charges 
couples to the local magnetic moment of the magnetic 
host, leading to the formation of a noncollinear magnetic 
state. This situation may be realized in lightly doped 
cuprate superconductorsi 4 i 5 i 6 i 7 i 8 i 9 i 10 ' 1:L i 12 

The main characteristic of the noncollinear state is 
that the spin configuration must be described by a set 
of three orthonormal vectors or, alternatively, by a ro- 
tational matrix which defines the orientation of this set 
with respect to some fixed reference frame. As a con- 
sequence, the order-parameter space is isomorphic to 
the three-dimensional rotational group 5*0(3), and in 
the low-temperature phase, when the rotational symme- 
try is fully broken, three spin-wave modes are present 



in the system, instead of two, as in the nonfrustrated 
case. Moreover, topological defects may arise in the sys- 
tem, associated with the chiral degeneracy of the spiral, 
which can rotate clockwise or counterclockwise. Because 
the order parameter space has a nontrivial first homo- 
topy group, 7Ti[50(3)] = Z 2 , the topological excitations 
are vortex-like. On the other hand, skyrmions are not 
present because the second homotopy group of the 50(3) 
is trivial, tt 2 [50(3)] = OM. 

A convenient field-theoretical description of frus- 
trated Heisenberg systems in the long-wavelength limit 
is provided by the 50(3) nonlinear sigma (NLc) 
modeliiiiiiiSiiSii 7 . Its critical behavior in two dimensions 
has been extensively investigated, both in the absence 
and in the presence of topological excitations. Studies in 
the former case have revealed a dynamical enhancement 
of the symmetry from 0(3) ® 0(2) to 0(4) under renor- 
malization group flow in d = 2 + e, which means that in 
the critical region all the three spin- wave modes have the 
same velocityiMiiZ. When topological excitations are in- 
cluded, a complex finite-temperature behavior is founds 
Numerical studies, as well as analysis involving entropy 
and free energy arguments, indicate the occurrence of 
a transition driven by vortex-antivortex pairs unbinding 
at a finite temperature ,I2i22i2ii22i In contrast to the 
XY case, here vortices and spin-waves are coupled al- 
ready in the harmonic approximation, and anharmonic 
spin- wave interactions yield a finite correlation length for 
arbitrarily low temperaturesiii 2 ^ Therefore, the transi- 
tion mediated by vortices is rather a crossover than a 
true Kosterlitz-Thouless (KT) transition. 24 Free vortices 
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start to proliferate at the temperature T v , similarly as 
vortices in the XY model do above the KT-transition 
temperature. 

In the present paper we study the physical properties 
of frustrated Heisenberg spin system, which are sensi- 
tive to the dynamics of the above-mentioned topological 
defects. The approach we use has been employed to de- 
scribe the dynamics of excitations in a very broad class 
of one- or two-dimensional systems^ The central idea is 
the application of the collective-coordinate method^ to 
quantize a non-trivial static solution of the classical equa- 
tion of motion of the field-theoretical model in question. 
In our case, we find that single vortex-like excitations or 
vortex-antivortex pairs are the localized static solution of 
the 50(3) NLct model. A proper description of the quan- 
tum levels associated with these solutions is provided, on 
the semi-classical level, by a theory in which the topo- 
logical excitation is represented by a single quantum me- 
chanical variable coupled to a bath of quantum harmonic 
oscillators, which are the fluctuations about the classical 
solution itself. Thus, the resulting effective model repre- 
sents a particle (the topological defect) scattered by the 
linearized excitations of the system. The latter can be 
integrated out using the standard system-plus-reservoir 
approach^ leading to a dissipative equation of motion 
for the topological excitation. As a consequence, any 
physical property of these systems that depends on the 
motion of the topological excitations may be expressed 
in terms of transport coefficients - such as mobility and 
diffusion - of these damped defects. Since we do neglect 
any interaction between the defects, our results are only 
valid for a diluted gas of topological excitations. Part of 
our results concerning the mobility of a vortex-antivortex 
defect has been recently published in Ref. [T3- Here, be- 
sides a complete presentation of the technical details, 
we discuss also the transport at finite frequencies and 
we compare -qualitatively and quantitatively- the cases 
where the defect is represented by a single vortex or by 
a vortex-antivortex pair. Moreover, in the light of new 
experimental results by Ando et ai.jSl we discuss the rel- 
evance of the single vortices for the transport in the spin- 
glass phase of cuprates. 

The structure of the paper is the following. Starting 
from the 50(3) NLcr model, we derive in Sec. II the quan- 
tum Hamiltonian describing the dynamics of the topo- 
logical defect coupled to a bath of magnetic excitations. 
In Sec. Ill the equation governing the evolution of the 
reduced density matrix for the topological defect is ob- 
tained and the influence functional, which describes the 
effect of the magnon bath on the dynamics of the vortices 
is evaluated. Section IV is devoted to the derivation of 
the effective action for the defect after the magnons have 
been integrated out, and in Sec. V the inverse mobility is 
calculated. In Sec. VI we discuss how a spiral state may 
be realized in cuprates and we then apply our results to 
this specific case. Section VII contains our conclusions. 
Details of the calculations are given in the Appendices. 



II. THE MODEL 

In the spiral state the spin configuration S at each site 
r is described by means of a dreibein order parameter 
n a k G 50(3) with k = 1, 2, 3 and n a k n a q = 5 kq ^ so that 

g 

— = n = ni cos(k s ■ r) - n 2 sin(k s • r), (1) 

where 5 = |S| and the wave vector k s = (jr/a^/a) + 
Q, with a denoting the lattice constant. Here Q — 
(2n/m x a, 2ir/m v a) measures the incommensurate spin 
correlations . Indeed, the magnetic susceptibility corre- 
sponding to the spin modulation has two peaks at k s 
and — k s (equivalent to (w/a,iv/a) — Q), as represented 
in Fig. 1 in the case of m x = —m y . The resulting spin 
order for ni and n 2 in the plane is represented in Fig. 2, 
where m x = —m y = 20. Observe that the periodicity of 
the spin texture is 27r/Q for even values of m x , m y , and 
twice it for odd values. 
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FIG. 1: Incommensurate magnetic response for the spiral spin 
modulation Q . The magnetic susceptibility corresponding to 
the spiral order with the wave vector k 3 exhibits two peaks 
at the points (n/a, n/a) ± Q marked by a circle in the figure. 
In this case Q has finite components in both the x and y 
directions, and the distance between the peaks is twice the 
modulus of Q. 

As discussed in Refs. a proper continuum field 

theory for the spiral state is provided by the 50(3) quan- 
tum NLct model, 

S = J dtd 2 x[K k (d t n k ) 2 - p ka (d a n k ) 2 }. 

Here, the index a stands for the spatial coordinates and 
summation over repeated indices is understood. The 
spatial anisotropy of the spin stiffness p ka depends on 
the components Q a of the incommensurate wave vector. 
Since at the fixed point all the spin-wave velocities are 
equalji^ii we will consider the case K k = k, p ka = p a 
and we will choose a system of coordinates parallel (x||) 
and perpendicular (x±) to the spiral axis, respectively, 

S = k J dtdx ± dx\\[{d t n k f - cl(d ± n k ) 2 - cf (d||n fe ) 2 ], 

(2) 
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FIG. 2: Spin background corresponding to Eq. Q and 
to the case depicted in Fig. 1. Here a value of Q = 
(27r/20a, — 27r/20a) has been choosen. 



and the corresponding Lagrangian obtained from the ac- 
tion (JSJ reads 



U = jAf I d 2 r(d^a) 2 



(7) 



with d^Ad^B = d t Ad t B - c 2 VAWB. By making the 
ansatai^ a(r,t) — mty(r,t), where to is a constant unit 
vector and \& a scalar function, the Lagrangian (7J re- 
duces to 
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L = -TV / d 2 v(d^) 



The equation of motion for the field VE", 
d 2 ^ - c 2 V 2 * = 0, 



(8) 



(9) 



possesses static topologically non-trivial solutions in the 
form of a single- vortex defect at R = (X, Y), 



where c±_ = \/pj_/ n and cii = y/pu]~K are the spin wave 
velocities perpendicular and parallel to the spiral axis. 
Even though for the moment we will keep our derivation 
on general grounds, in Sec. VI we will specify the values 
of the parameters k and c for the case of cuprates, where 
they can be related to measurable quantities. 

Given the action J5J as our starting model, our first 
aim is to analyze whether the equations of motion admit 
topologically non-trivial solutions. For that purpose, it 
is convenient to introduce an equivalent representation of 
the rtf, order parameter through an element g e SU(2) 
as n£ = (l/2)tr[a a ga k g~ 1 ], where a a are Pauli matrices, 
and to introduce the fields 



4" 



1 



2i 



(3) 



which are related to the first derivatives of the order pa- 
rameter through d^nl — 2a jkA^n"-^ Here, = (d t , V) 
and £i23 = e 123 = 1. Using that (<9 M n) 2 = 8A 2 (no sum- 
mation over index [i is imposed here) , the action J2J) reads 



S = 8k / dtdx\\dx\ 



4Ai 



„2 A 2 



(4) 



The above action may be mapped to an isotropic form 
by introducing the coordinates r = (x, y) with 



V = 



We then find 



S = M J dtd 2 rA 2 t = M J dtd 2 r (A 2 - c 2 A 2 ) 



(5) 



with the isotropic spin-wave velocity c = ^/C||Cj_ and the 
constant TV" = 8k. The most generic expression for the 
element g is given by 



g[a] = exp ( ~a(r, t) ■ a 



(6) 



'Fit, = arctan 



and a vortex-antivortex pair 



x-X 
y-Yj ' 



(10) 



^2v = arctan 



= arctan 



— arctan 



x-Xx 
y-Yi 

[dx(r-R)], 
[(r-R) 2 -d 2 /4] 



x-X 2 
y-Y, 



(11) 



where now R = (Ri + R2)/2 is the center of mass and 
d = R2 — Ri the relative coordinate of the defect pair. 
If in Eqs. (|10|) and the role of x and y coordinates 
is interchanged, one only changes the vorticity. Without 
loss of generality, we may assume the unity vector to be 
in the z-direction, m = e z . Thus, the rife fields which de- 
fine the spin configuration according to Eq. <£Q), are given 
by m = (cos*,-sin*,0),n 2 = (sin 'F, cos >F, 0), n 3 = 
(0, 0, 1). The spin patterns corresponding to a single vor- 
tex (<F from Eq. 1)100) or to a vortex-antivortex pair 
from Eq. are represented in Fig. 3 and Fig. 4, re- 

spectively. 

The main difference between the two possible static so- 
lutions l|10|) and i|ll|) is their energy. As shown in App. A, 
the energy of a single- vortex diverges with the logarithm 
of the system size £, E[*i>i v ] cx ln^?. On the other hand, 
the vortex-antivortex pairs have finite energy, depending 
on the distance d between defects, -E[5 , 2d] oc lnc?. A simi- 
lar situation is realized in the standard XY model, where 
indeed the presence of single defects below the KT transi- 
tion is not energetically favorable in the thermodynamic 
limitiSi However, in the case of the model (J2J, which 
possesses asymptotic freedom, the correlation length £ 
is finite at any finite temperaturefiiS so that the loga- 
rithmic divergence of the single-vortex energy should be 
understood up to the length scale £, £[^1,,] oc ln£. In 
addition, the energy of the vortex-antivortex pair should 
be bounded below at distances of the order of few lat- 
tice spacings, which is the intrinsic cutoff of the theory. 
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FIG. 3: Spin background corresponding to the one-vortex 
solution of Eq. HUH . The center of the vortex is marked by 
the circle. The spiral incommensurability Q is the same used 
in Fig. 2. 
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FIG. 4: Spin background corresponding to the vortex- 
antivortex solution of Eq. 111H . The centers of the vortices 
are marked by a circle. The spiral incommensurability Q is 
the same used in Fig. 2. 



Because the procedure that we describe in the following 
does not depend on the exact form of the static solution, 
we will refer to a static topological defect solution ty v 
specifying only at the end of the calculations the differ- 
ences between the cases (|10|l and (lllll . 

Follo wing a procedure analogous to the one described 
in Ref. |26| to quantize the kink solution of the scalar 
field theory, we analyze now the effect of the fluctuations 
around the static topologically nontrivial configuration, 
which is a saddle point of the action that corresponds to 
the Lagrangian (JSJ). In order to reach this aim, we write 
the generic field g € SU (2) of Eq. © in the form of a 
product of the field g s corresponding to a static solution 
rn^> v (r) and the field g £ corresponding to the fluctuations 
around it, 



where 



g e [e\ = exp ( -e ■ a 



Observe that the description of the fluctuations via Eq. 
(|12|l differs from the standard approach used for a scalar 
field theoryj2& and it is related to the symmetry prop- 
erties of the order parameter. Indeed, since the full 
g has to be an element of the SU(2) group, and both 
g Sl g S SU(2), then the fluctuations g s around g s have 
to belong to SU{2) as well. If instead, we had used the 
expansion a = rn 1 i' v + e, the equations of motion for the e 
field would have been independent of the static solution 
leading to a failure of the semiclassical expansion. 
Using Eqs. © and (|T2*)l . we can express the action © in 
terms of the fields ^ v and e. Retaining only terms up to 
second order in e, we find (Appendix B) 



1 / e 2 + e 2 



e a s b m b d^ v 



+ \e abc e b m c d^ v + ^e a + ±e abc e b d^e c , 

where e 2 = e 2 + e 2 . The corresponding Lagrangian then 
reads (App. B) 



L +N 



d 2 rd 



(13) 



with L given by Eq. JSJl and 

A[*i>] = \(d^e) 2 + \{m-d ll e)d^ v 



- A e abc d tl e a e b m c d^ v 
1 



-(d,e) 2 + -e 2 (d ll 6) 2 + -(d,e z ) 2 



~£ 2 a/i>,.,A 



(14) 



Here, we used the fact that m = e z and introduced polar 
coordinates e = (ecos9,esm8,e z ). Since the Lagrangian 
L\ is evaluated at the vortex- like solution ^ v of Eq. 
the equations of motion for the fluctuations around the 
topological defect also depend on ^f v 

e : (d 2 - c 2 V 2 )e - e{d^6) 2 - c 2 eW„V0 = 0, (15) 



dJs 2 d u 



-V(e 2 W„) = 0, (16) 



(d 2 - c 2 V 2 )e z - 



(17) 



Eq. ijToj) admits the solution — ^ v /2, whereas Eq. l(T7|) 
indicates that the field e z is free. By using these two 
conditions, we can rewrite the total Lagrangian L in Eq. 
CCU as 



g(r,t)=g a [tf„(r)] g £ [e(r, t)} , 



(12) 



{d^ v ) 2 + {d,s) 2 --e 2 {d^ v f 



(18) 
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and the equation of motion (II 5|) as 

1 



e = 0. 



where the last sum is over all non-zero-frequency modes. 
By inserting these expressions into the full Lagrangian 
(TL%1) evaluated at the saddle-point solution we obtain 



Since the field ^ v in the previous equation does not de- 
pend on time, we decompose the field e into its time- and 
space-dependent parts 



Y / dtd?r(d t * v ) 2 = — I ,ll<l-r II, R, <>„*,J>M, 



A I 



dm 2 (t), 



(25) 



(19) 



where 



and identify the normal modes rj nm with the eigenfunc- 
tions of the operator 



M 



d 2 r(W„) : 



c 2 [V 2 + V(r)] J] nm = -w 2 ™^ 



(20) 



This equation has the typical form of a Schrodingcr-likc 
equation for a particle scattered by a potential V(r) — 
(V^ v ) 2 1 '4. The two indices n and m refer, respectively, 
to the radial and angular part of the wavefunction. By 
using a standard approach to scattering problems in two 
dimensions one may express the wavefunctions ry nm in 
terms of the eigenfunctions of the free problem (V = 0), 
corrected by a phase shift 5 m due to the scattering by 
the potential V(r)^ 



is the mass of the topological defect, which is propor- 
tional to its energy (see Appendix A). The time deriva- 
tive of the field e yields 

^ J dtd 2 r(d t e) 2 = X E / dtd 2 r[q* nm (tW nm 

nm,kl 

- Q*n m d a V* nm Ra{t) + C.c.j [q k i(t)r] kl - q kl d[3T] kl R0(t) 



+ C.c. 



kl 



1 



2f 



H^\(k nm r)+e- 2 ^HlZ(k, 



r(2), 



| m | 



I m I 



(1,2) 
|m| 



(21) 

are Hankel functions of the first and sec- 



M ) 



.kl 



(26) 



Here, 

ond kinds, m is an integer, and $ is a polar angle. The 
k nm values are determined by requiring the vanishing of 
the wavefunction l|21|l at the boundary r — £. By using 
the asymptotic form of the Hankel functions we obtain 
knm£ = (2n + l)7r/2 + (2|m| + 1)tt/4 — S m , where n is a 
positive integer. Since the field e is real, we may rewrite 
the expansion l|19|) in the form 



where the coupling constants G are related to the eigen- 
functions rj via 



M 



d 2 rr ]kl Vr 1 * n 



and we neglected terms of order g 2 R 2 . Here, we used 
that / d 2 vr\ nm r\ k i = and / d 2 rrj nm Vi] k i = for m and 
I positive. By substituting Eqs. Ij25(l and Q26|l into the 
Lagrangian 118|) . we obtain 



e(r,*)= J2 (lnm(t)Vn m (r) + q* m {t)Vn m ^)) , (22) L = hlK 2 + ^(|g„ m | 2 - u; 2 nm \q nm \ 2 ) 



n,m>0 



where we used the identities rj n>m = e -2 m ry*_ m and 
$m = Note that the sum in Eq. 12 21) is over the 

positive angular momenta, as will be the case in what 
follows. 

The static defect solution ^ v (r) of Eq. 10 is invariant 
under translation of the center of the defect (i.e., the po- 
sition of the vortex or the center of mass of the vortex- 
antivortex pair). A consequence of this invariance2& is 
that Eq. I|2U|) admits zero- frequency modes. A consistent 
treatment of them requires the use of the collective co- 
ordinate methodi2Si2& The center of mass of the defect is 
then promoted to a dynamical variable, yielding 



7~1 R(*) (9nm9«G* m kl + q* m q k iG 



nm,kl I ) 



nm,kl 



where we rescaled q — ► q^jN/2. Using that G* m kl = 
— Gkl,nm an d neglecting terms which are quadratic in G, 
the corresponding Hamiltonian reads 

H = -i-(P - V E ) 2 + ^(| Pnm | 2 + ul m \ qnm \ 2 ), (27) 



with 



* (r)->*„(r-R(i)), 



(23) 



and 



e(r,t) -> e(r-R(f),t) = ^( ?nro (t)?7„ ro (r-R(t))+C.c.), 

nm 

(24) 



nm,kl 

Here, P is the momentum canonically conjugate to the 
center of the defect R, and q nrn and p nm are the co- 
ordinates and momenta of the magnons. The classical 
Hamiltonian l|27|) can be promptly quantized by intro- 
ducing two sets of independent creation and annihilation 
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operators, a), a, and 
reads 



b. The quantum Hamiltonian 



H — H v + Hb + Hi , 



where 



H v 



2M 



is the Hamiltonian of a free defect, and 

ii^nm^nm + b nm bnm) 



(28) 



(29) 



(30) 



is the Hamiltonian of the bath of magnons which consists 
of two independent sets of noninteracting harmonic os- 
cillators described by the operators a, and b, b', as it 
is expected in two dimensions. The interaction between 
the bath and the topological defect is described by the 
Hamiltonian 



Hi = -j-jr ^ ^ |^D nm ^^Ci nm fl^ ^kl,nrnbnmb\i 
nm,kl 



aklbnm) 



(31) 



with the coupling constants given by 



^nm,kl 
^rt7n,kl 



2 y/iJ. 

1 L^nm 



ara^kl 
Ukl 



2 y/UJ 



Qnm.kl ■ 



(32) 



The terms with the coupling constants D/C describe 
the scattering/creation (annihilation) of magnetic excita- 
tions by the defect. Since we consider here only the low- 
energy dynamics of the topological defect, we neglect the 
off-diagonal terms in the interaction Hamiltonian, i.e., 
we set C = 0. In the following we shall integrate out the 
bath degrees of freedom in order to study the effective 
dynamics of the defects. For that purpose we employ the 
Feynman- Vernon formalism. 



III. REDUCED DENSITY MATRIX 

In this section we derive the reduced density matrix 
for the defect. The system under consideration consists 
of two subsystems: the topological defect and the bath 
of magnons. Thus, the Hilbert space of the full sys- 
tem, H, is a direct product of the subsystem Hilbert 



spaces H — H v <g> Hb = H v <g> H 



(«) 



Hg\ and the 
state of the full system is also a direct product, |x, d?) = 
|x) eg) \a) = |x) ® \a a ) (£> |d?b). We use the coordinate 
representation for the defect (x are the eigenvalues of its 
center-of-mass position operator) , and the coherent state 



representation for the bath, a nm \d a ) = a., 



and 



b n m\db) — ctnm,b\db)- The reduced density matrix is de- 
fined as p v (t) — trs[/5(i)], where trg denotes the trace 



over the bath degrees of freedom, and p(t) is the density 
matrix of the full system, whose evolution is described 
by 



iHt iHt 

pit) = e * p(0)e h . 



(33) 



Here, H is the Hamiltonian of the full system given 
by Eqs. (|28|l - l|32|) . The matrix elements of the density 
operator in the basis introduced before are 

P(x,a;y,/3;i) = (x, a\p(t)\y,/3), 
and the reduced density matrix of the vortex reads 

d 2 d 



d 2 a , 



T 2N 



(x,a|/5(t)|y,a) 



T 2N 



x,a|e 1 " t / 5(0)e 1 ^ i |y, a) 



(34) 



After insertion of the unity operator on both sides of p(0) 
in Eq. H34[) , the reduced density matrix acquires the form 



7T 

■Ht. 



a 

2JV 



d 2 x' / d 2 y' 



(x, a\e-^ |x', /3)(x', /3|p(0)|y', f3')(y', f3'\e^ |y, d). 

(35) 

In order to calculate the time evolution of the reduced 
density matrix, we have to define the initial condition for 
the density matrix of the full system. For the sake of 
simplicity, we choose the factorizable one 



p(0) = K(0)Pb(0), 



(36) 



which implies that the bath and the topological defect are 
decoupled at t — 0. The bath is assumed to be initially 
in thermal equilibrium at temperature T, 



Pb 



,-UH b 



tr[ e -UH B ] 



— PB,aPB,bi 



(37) 



where U = ft/(fcsT). Here, we used the fact that the 
baths do not interact, so the density matrix of the full 
bath is the product of the density matrices for the sepa- 
rate baths. By substituting Eqs. i|36|) and l|37|) into Eq. 
(|35|l we obtain 

X(x,y,t)= J rf 2 x' J d 2 y'J(x,y,i;x',y',0)^(x',y',0), 

(38) 

with the superpropagator J given by 

J(x,y,t;x',y',0)=y — J -J^ J -J^p^* , ? , 0) 
A-(x, a* ; x', 0; t)K* (y, at ; y', 0'; t). (39) 



J 
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A. The Superpropagator 

We consider the superpropagator (|3"5j) . The kernel 

K(x,a*;y,0;t) = (x, ^e" 1 * 1 |y, 0), 
can be expressed in the path-integral formalism as (Appendix C) 

12 \(3\ 2 ) r*„- n 



if(x, a*; y, /3; t) — / Pxexp 



|a| 
2 2 



PCexp <j -[C*(0) • + C(t) ■ «1 > exp j -S [x] } exp 



where 5o[x] = / * dt'(M /2)ic 2 stands for the free action and 

s/[x, cl = jf* { -of • - c ■ t) - \ [mc\ o x . hj(c : c) 

denotes the interaction. Here 



Sl, a [x, Ca] + <S/,&[x, Cfc] 



(40) 

{^/[x,C]}, 
(41) 

(42) 



— J]] Fl^nmCnm^Cnmd — ^Bm + hs,b, hi — H ^nm,klCkLaCnm,a ~ 'Dkl,nmCnm,bCkLb = h-J.a + h/,i 

nm,i— a,fe 



nm,kl 



By inserting Eq. (|41|l into Eq. 139|) and using the reduced density matrix of the bath in the coherent state representation 
(Appendix D), we obtain a simpler expression for the superpropagator, 



J(x,y,i;x',y',0) = fpx f Vye^ s ° ^~ s «^F{x, y], 



(43) 



where = J- a Tb is the total influence functional, and Ti (i — a, b) is the influence functional for the bath i given by 
(to simplify notation we omit the index i in the integration variables) 



d 2 a f d 2 f d 2 0' % 



VQ / Pfexp - C*(0) • + C(t) ■ a* + 7(0) ■ 0'* + ■ « 

/3 7/3'* I * 



exp 



{^:[x,C]+%[y,7]}, (44) 



with the initial conditions 



aO)=0, C(t)=a*, 



(45) 



T(O)=0'* i(t) = a. 



(46) 



B. Influence Functional 



We now evaluate the influence functional, which describes the influence of the bath on the effective dynamics of 
the defect. The only difference between the functionals T a and Tb is in the form of the interaction Si; see Eqs. (|42f) 
and (|44|l . Note that the actions S7,a and S7,!> are related by the substitution ~D nm ,kl ~ * — ^H,nm- Thus, it is enough 
to calculate the functional T a , and consequently T\, is obtained using the latter transformation. In order to simplify 
notation, in what follows we write the integration variables without the index a. First, we calculate the path-integrals 
in Eq. I|44|) using the stationary phase approximation (SPA). In order to apply the SPA, we have to solve the equations 
of motion corresponding to Si and Sj. Because Si(x, () = — S/(x, £), we need to consider only Si. The equations of 
motion are promptly obtained from 5Si/5(* m — 0, SSi/S£ nm = 0, and they read 

Cnm "f" iLJnmCnm ^X ^ "^kl,nmCkl 0; 

kl 

Cnm ~ i^nmCnm + *X ~^nm,klCkl = 0" 
kl 

(47) 
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Notice that the two equations are identical, one is the complex conjugate of the other (recall that ~D nm .ki — ). 

The SPA requires the evaluation of the action Sj on the classical trajectory, which is the solution of the above 
equations of motion. Straightforward calculations show that the value of £/ at the stationary point is zero. If we 
define £ = ( c i + A, then the functional integral over £ becomes the functional integral over the fluctuations A around 
the saddle point. Expanding the action around its saddle point, we find that the relevant contribution comes from the 
second derivative of the action at the stationary point, because both the value and the first derivative of the action are 
zero at the saddle point. The second derivative of the action evaluated at the stationary point is a constant operator, 
so the integration over the fluctuations A yields, up to an irrelevant constant, 



d 2 a 



F a [x,y] = I Z£ I ±P J ^ pBA ^jr )e 



exp<;- 



C*(0)-/? + C(t)-a*+7(0)-/3'*+7*(t)-« 



(48) 



Therefore, in the SPA, Si and Sj only contribute to the 
influence functional through the boundary terms, which 
may be determined using the solutions of the equations 
of motion (see Appendix E). 

After inserting Eqs. (|D2|) and (|E8|I into Eq. (gHJ), and 
performing the Gaussian integrals over a, (3 and /3', the 
influence functional acquires the form 



Fa [x. y] 



det(l -nT a Y 
where the matrix T a is given by 



(49) 



r£ ro , w [x,y] = -[W H , nm (x,t) + WWfcz(x,0) 
+ W u , nm (y , 0) + W nmM (y,t)] + _^ [W nm ,ki (x, 0) 

pq 

+ W pq , nm (x,t)}[Wki,p q (y,0) + W pqM {y,t)], (50) 

and n pq = l/[exp(Uuj pq ) — 1] is the bosonic occupation 
number. Equation l|E7|l enables us to express the matrix 
r a only in terms of the functionals W and W 

r" m ,fei[x,y] = i[W fcZ , nm (x,t) + W nm ,idfc°) 

+ Wu,nm(y> °) + W: mM (y, t)] + 1 ^[f„ m , M (x, 0) 

pq 

+ W pq , nrn (^t)][W^ pq (y,0) + W; qM (y,t)]. (51) 
Using the formula lndet„4 = trln.4 for the matrix A — 



(l-nT ) -1 , we find 



•F a [x,y] = exp[tr(nr Q )] = exp 



pq,pq 



(52) 



and the total influence functional reads 



T = T a Tb = exp 



in the lowest order in nT, 



L pq 
where T 



l>pq^ pq,pq 



(53) 



r Q + r & . The diago- 
nal elements of the matrix T a are obtained from Eq. Q51f). 
while the matrix T b is obtained from T a by the substi- 
tution T> nm ,kl — > -T>kl, nm = -^nm,kV Tlie f unctionals 
W and W are given implicitly by Eqs. (|E3|) . From their 
form we see that they actually represent the amplitude 
of scattering of the mode nm to the mode kl through 
virtual intermediate states. These functionals can be de- 
termined iteratively from these equations up to any or- 
der. Here, we study the motion of a vortex with small 
kinetic energy; therefore the Born approximation will be 
enough for our purpose. The functionals W and W are 
calculated within the Born approximation in Appendix 
E. Using Eq. i|E4|) the diagonal elements of the matrix T 
can be promptly evaluated (Appendix F), and the total 
influence functional reads 



^[x,y] 



exp[-$] exp[<£>] 



(54) 



where 



2 pt pt 

Yl / dt l dt"8(t' -t")e» v (t' -t")[x»{t')-y»(t')][x v (t")+y 1 ' {t'% (55) 



fJ>,U— 1 



JO 



$ = d£ [ dt''8{t' -t'')^%t' ~t l ')[x^{t l )-y^{t')][± u {t l ')~y v {t 1 ' 

H» JO JO 
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with 

e""(t) = E (C,^«,« + £C jW ^,«) M"nm ~ w w )t, (56) 

nm,kl 



nrn.kl 



From Eqs. I|43[) and l|54|l we see that the oscillatory part 
exp[i$/fi,] gives a contribution to the effective action of 
the defect due to its scattering by the magnons and leads 
to its dissipative motion, as we show in the following sec- 
tion. The decaying part exp[$] is related to the diffu- 
sive properties of the vortex. The diffusive and damping 
properties of the defect are related at low temperatures 
by the fluctuation-dissipation theorem. 



IV. THE DYNAMICS OF THE DEFECT 

A. Transport properties of the defect 

In this section we shall study the effective dynamics of 
the defect after integrating out the magnons. According 
to Eqs. i|43|) and l|54l) the effective action describing the 
influence of magnons on the motion of the topological 



defect reads 

S , o// = 5o[x]-Sb[y] + *[x,y], (57) 

where $ is given by Eq. (|55|) . Since $ a D we observe 
that if the coupling constants D were zero, then the mo- 
tion of the defect would be free. The equations of motion 
for the defect can be directly obtained by extremizing the 
effective action l|57|). 5S ef f/5x^ = and SS e f f /6y^ = 0. 
In terms of the center of mass v = (x + y)/2 and relative 
u = x y coordinates, they read 

^(t)+E/ dt'^ u {r-t')v v {t')=Q, 

^(t)+E/ dt'^ v (t' -T)ii u (t') = 0. (58) 

As a consequence, the damping matrix r y f11 ' is given by 



nrn.pl 

~A~h7 2-^ [ - nnm ~ n P l > 7 7 \ Lr nm,pl Lr nm,pl + ^ nm ,pl U nm ,pl ) COS(U> nm - U! p l)t 

*±1V1 UJ nrn UJpl \ / 



nra.pl 

= V Km ~ n pl) " — G nm pl G run,pl COs(w„ m - LJ pl )t, (59) 

ZlVl — ' UJ nm U! p i 
nrn^pl 

where e^ v was defined in Eq. J5BJ. Introducing the scattering matrix 

S"> jW ') = E - W»m)%' " M (60) 

tim.p! 

and observing that because of the isotropy of the model the damping matrix is diagonal (see also Eq. H70|) below), 
•yV(t) — 7(t)<5 A " / , the damping function j(t) is given by 

7 (t) = _JL fd, /°° <fc/[n(tj) - ft(h/)] ^ ~ + ^ g(u>, u/) cos(cj - u% (61) 



Let us introduce the new variables £ = (w + w')/2, ip 
lj — lu' . The damping function can be rewritten as 



where J~(ip) is the spectral function of the bathS given 
by an additional integration 



7(<) 



dipj{ip) cos V>t, 



<62) - -sr* 
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£ 2 - V 2 /4 



(63) 



From the equations of motion (|58|l it is easy to see that if 
a charge q is associated with the defect (see next section) 
the corresponding optical conductivity is 



nq 2 /M 



■7(2) 



(64) 



where n is the density of carriers and 7(2:) is the Laplace 
transformation of the damping function 



%z)= / dt 7 (t)e- 



dip 



z 2 + %p 2 



(65) 



When only quasielastic processes are taken into account 
ip w 0, so that we can approximate i7(V0 ~ »7(0), and 
using the fact the dtpcosipt = 2nS(t), we find from 
Eqs. (E2J> and JSHJl that 

7(<)=7<*(*). 7(^=7/2, 
with the damping coefficient 7(T) 



7(T) = 2nJ(0) 
where we defined 



47T?l 



gn(g) 
9^ : 



(66) 



(67) 



According to Eq. I|64|) the real part of the optical con- 
ductivity has a typical Drude-like shape 



Rea{<jj) 



nq 



7/2 



M oj 2 



(7/2) 2 



(68) 



where 7/2 plays the role of the inverse scattering time. 
It is then natural to introduce the inverse mobility as 



M7 
~2q~- 



(69) 



It is worth noting that even though the formula 164(1 
is general, in the computation of the damping func- 
tion (|59() , we considered only low-energy quasielastic pro- 
cesses, which, naturally, should lead to the Drude-like re- 
sponse (|68p. If one keeps in the evaluation of J{ip) the 
next nonvanishing contribution to Eq. I|65l) . 

Jty) «7(0) + J 2 1p 2 , 1p<LO Cl 

where lu c is a proper cutoff for the previous expansion, 
7(z) can be estimated as 

J(Z) = 7T [,7(0) - z 2 J 2 ] + 2J 2 LO c Z. 

The corresponding optical conductivity is then 
nq 2 wj(0) + ttJ 2 u 2 



Rea(u>) 



M w 2 (l + 2J 2 oj c ) 2 + (nj(0) + itJ 2 uj 2 ) 2 ' 



which is qualitatively the same as the one given in Eq. 
B68JI . In particular, the dc conductivity is found to be 
the same. Since in the following we will address the issue 
of the temperature dependence of the resistivity, we can 
safely rely on the approximation l|66l) of the damping co- 
efficient, which takes into account only the contribution 
of ,7(0). 



B. Evaluation of the damping coefficient 

In order to determine the damping coefficient l|6T?|) , we 
first have to evaluate the function <p(£) defined in Eq. 
()67J) . By rewriting the summations over the (radial) in- 
dexes n,p in Eq. Ij60|l as integration over continuum vari- 
ables, and taking into account that Aq n ~ n/£ we find 
that 

S{k,k') = S yy {k,k') = S xx (k,k') 
= ^Hj d( lj dq'\G x qmtq , l \ 2 6(k~q)S(k'-q / ) 



ml 



1 2 

m,k'l\ 5 



nil 



S xy (k,k') 



}^tm.k'l^km.k'l — 0' (70) 



where k — to/c, k' = lo'/c and the last equation follows 
from Eq. i|G8|) . The limit l|67|l thus acquires the form 



p(0 - ~ 2 /h Urn (k k') 2 \G x mk ,lv f- ( 71 ) 

Till 

From the above relation it is obvious that the only terms 
of fe , ; | 2 which contribute to the limit l|67|l are those 
behaving like ~ l/(k — k') 2 , i.e., the term proportional 
to A^ 2 ) in Eq. (JSEJ. We then find 



2fc 2 



Vir (2) 1 2 = 

Z^,' km,k'l\ P(k — k') 2 

ml y ' 



where 



A 



(2) 

m+l,m 



' e «/2 + e -i7r/2 e 2i(a m+1 -5 m )^ / e -i7r 2 



e i7r/2 e -2i(a ro+ i-5 m ) 



4 sin 2 (6. 



m+1 "mi 



(73) 



Substituting Eq. (|73|l into Eq. (|72|l . and evaluating the 
limit H71f) we finally obtain 



¥>(0 = 4^e(0, 

where £ = cfc and we defined 

00 

5 = E sin 2 (5 m+ i - (5 m ). 

m=0 



(74) 



(75) 
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Using the preceding equation, the damping coefficient 
acquires the form 



l{T) 



32h 2 



Mnc 2 kBT J 



^ e k B T 



( h e 
I e k B T 



" ^ \ 2 ' 



(76) 



Observe that Eq. (|76|) is valid for both kinds of defect 
solution (|10fl and lllj l. However, since the phase shifts 
are determined by the eigenfunction rj of the scattering 
problem 120(1 . or will give rise to different po- 
tentials Vi„(r) or V2n(r), and then different phase shifts. 
As a consequence, also the function G{Q in Eq. 1(75(1 will 
be different in the two cases, leading to a different tem- 
perature dependence of the damping coefficient ((76(1 . 



V. INVERSE MOBILITY 

We evaluate the phase shifts by adopting the Born 
approximationi^&Si The phase shift S m (k) of the wave- 
function with angular momentum m and wave vector k 
then reads 



S m (k) = arctan7r^4 m (fc), 



(77) 



where A m (k) is the expectation value of the potential 
over the eigenfunction of the corresponding unperturbed 
Schrodinger equation, i.e., the Bessel function J m (kr) in 
the case of Eq. PU|) 

/■oc 

A m {k)= / drr[J m (kr)] 2 V(r). (78) 
Jo 



A. Vortex-antivortex pair 

Let us first consider the case of the scattering of a 
vortex-antivortex pair by the magnons. In this case the 
potential V(r) in Eq. 1(20(1 has the form 



V 2v (r) = i (V* 2 ,) 2 , 



(79) 



where ^ 2v = - with & lv given by Eq. (X2)| 

One can readily show that 



(r-Ri, 2 ) 2 ' 
(r - Ri) ■ (r - R 2 ) 
(r-R!)2(r-R 2 )2' 

which gives, using the translational invariance r — R — * r 

{ 2v) (r-d/2) 2 (r + d/2) 2 (r 2 + d 2 /4) 2 - (r • d) 2 

Since the distance d between defects is a fixed parameter 
and we consider length scales r > (i, we can approximate 
the potential as 



V 2v (r) 



d 2 /4 



(r 2 + d 2 /4) 2 



(80) 



For the potential ((80(1 , Eq. ((78(1 can be evaluated analyt- 
ically for m > 1, yielding 



A m (k) 



■nkd 



{I m (kd/2) [K m -i(kd/2) + K m+1 (kd/2)} 



K m {kd/2)[I m ^(kd/2) + I m+1 (kd/2)}} , 



where and K n are the modified Bessel functions of the 
first and the second kinds, respectively. 

Returning to Eq. (|7SJ), we observe that G(£, — k/c) 
is a function of the dimensionless variable y = £d/2c. 
Introducing the same variable also into Eq. 1(76(1 , we may 
rewrite the inverse mobility given by Eq. ((69(1 as 



_il28 E, 
Mo 



(81) 

where /i = h/ea 2 is the quantum of inverse mobility 
for a given lattice spacing a, d — aa and E c — hc/akB 
is the characteristic temperature scale associated with 
the magnons. Even though a quantitative estimate of 
Eq. 1(81(1 requires the knowledge of the values of these 
microscopic parameters, its qualitative behavior can be 
promptly understood. In particular, since all the tem- 
perature dependence of pT 1 is due to the Bose factor in 
Eq. ((81(1 . one can expect that the inverse mobility van- 
ishes at zero temperature, where no thermally activated 
scattering processes exist, and increases linearly at high 
temperatures, with the slope determined by the shape 
of the function G(y). In the left panel of Fig. 5 we plot 
(/x//xo) _1 a s a function of T/E c for several values of a. 
One observes that already at small fractions of the ra- 
tio T/E c , the inverse mobility is linear in temperature. 
Moreover, as a increases the linear behavior arises at 
even smaller temperatures, and the overall value of the 
inverse mobility decreases. 



B. Single vortex 

Let us analyze now the behavior of the inverse mobility 
obtained when we identify the defect as a single vortex. 
In this case, instead of Eq. l(79"|) the scattering potential 
reads 

As a consequence, the phase shifts defined by Eqs. I(77|) 
and 1(78(1 are given by (m > 1) 



5 m = arctan . 

8m 



(82) 



Note that the phase shifts in the case of a single- vortex do 
not depend on the wave-vector, but only on the angular 
momentum. Thus, the function G defined by Eq. 17511 . 
does not depend on the frequency 



64tt 2 



64tt 2 + [tt 2 + 64m(m+ l)p 



0.032. 



12 
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T/E c T/E c 

FIG. 5: (Color online) Inverse mobility in units of /i^" 1 as a 
function of the rescaled temperature T/E c . Left panel: in- 
verse mobility of a vortex-antivortex pair, according to Eq. 
I81H . at several values of a. Right panel: inverse mobility of 
a single- vortex defect, according to Eq. 183H 

Because Q is a constant, we can introduce the rescaled 
variable y = h^/ksT into the expression ilt)9|) for the 
inverse mobility, which yields 

^ {T ) = ^g[0 . (83) 

Here, we used the fact that /°° dyy 2 / (e y — l) 2 — tt 2 /3. In 
comparison with the case of a vortex-antivortex pair, the 
main difference is that here the inverse mobility depends 
on the square of the temperature, for all the tempera- 
tures. In the right panel of Fig. 5 we plot (nf/io)^ 1 as 
a function of T/E c : notice that the overall variation of 
the inverse mobility is smaller compared to the case of 
the vortex-antivortex pair, but they are still of the same 
order of magnitude. 

VI. APPLICATION OF THE MODEL: THE 
CASE OF CUPRATES 

A. The spiral state in cuprates 

Although the model we have developed above could be 
applied to describe the dynamics of topological defects in 
several frustrated Heisenberg spin systems, here we con- 
centrate on lightly doped cuprates. Indeed, a large part 
of the literature devoted to frustrated spin systems is con- 
nected to the t — J model, which is the strong-coupling 
limit of the Hubbard model. The latter is considered to 
be the prototype of an effective description of the CuC-2 
planes of cuprate superconductors. At half-filling the 
t — J model describes a spin- 1/2 antiferromagnet which 
is believed to have long-range order at zero temperature. 
As the system is doped away from half-filling the motion 



of a hole will leave a trail of spins pointing in the wrong 
direction. Thus, two issues must be settled: (i) the char- 
acter of the quasiparticle wave function and (ii) the effect 
of the hole motion on the spin background. Both issues 
have been extensively addressed in the literature. In the 
atomic limit (t = J = 0) of the Hubbard model, consid- 
ered by Brinkman and Rice^S the "string" of perturbed 
spins can be healed only by retracing the original path 
and returning all the spins to their original position. As 
a consequence, in the presence of a finite but small J, 
they argued that the ground state of the hole involves a 
magnetic polaron: the cost of creating a ferromagnetic 
region around the hole is compensate by the fact that in- 
side this ferromagnetic cloud it can sit at the free-particle 
band edge. 

At larger value of J (1 > J/t > 5 x 1CT 3 ) Shraiman 
and Siggia^ showed that at least in the Ising limit ( J z ^ 
0, Jj_ = 0) the picture of band narrowing effect is more 
appropriate than the polaron formation. In the Ising 
limit the holes are infinitely massive, because they are 
self-trapped to their original position by the "string" of 
overturned spins. When a finite J±_ is included, quan- 
tum spin fluctuations associated with it can repair a 
pair of overturned spins and the mass of the holes be- 
comes large but finite i^l Several calculations22i2&2i were 
performed using an effective Hamiltonian which cou- 
ples the holes (constrained to no double occupancy) to 
Holstein-Primakoff spin-waves. The result is that the 
hole moves on a given sublattice, forming a narrow (~ J) 
quasiparticle band with the minimum at the wave- vector 
(±7r/2, ±7t/2), plus an incoherent part originating from 
the spin-wave excitations created by the hole motion. By 
using a semiclassical approach, Shraiman and SiggiaA^ 
showed that one can assign to the hole states a dipo- 
lar momentum p a , which is a vector both in lattice and 
spin space. The coupling between this dipolar moment 
and the magnetization current j = SI x <9 Q f2 of the an- 
tiferromagnetic background, described by a NLrr model 
for $7, leads, at finite doping, to a spiral re-ordering 
of the antiferromagnetic phase of the background spins. 
Within a similar approach, Gooding^ has argued that 
a strong localization of the hole could eventually lead to 
a skyrmion-like configuration of the background spins. 
However, while the polaron or the skyrmion formation 
seem to be plausible scenarios for a single defect, at fi- 
nite doping the picture proposed by Shraiman and Siggia 
of a new helical spin configuration is more likely. Later, 
many calculations on the t — J (or t — t' — J) model based 
on different approaches have indeed confirmed that a spi- 
ral ground state can be favorable at low doping£2*2*i2i2ii 

Recently, the interest on spiral formation in t — J- 
based models has been revived due to the experimen- 
tal observation of incommensurate spin correlations in 
cuprates, i.e., an enhancement of the spin susceptibil- 
ity at a wave vector k s slightly displaced with respect 
to the commensurate wave vector (ir/a, n/a). In par- 
ticular, detailed measurements of the incommensurabil- 
ity as a function of doping are available for lanthanum- 
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based compounds For doping x larger than 0.05, 
i.e., in the regime where the samples are superconduct- 
ing, the observation of four peaks at (tr/a ± 5, 7r/a) and 
(•7r/a, it /a ± 5), and the simultaneous measurement of in- 
commensurate charge peaks, lead to a natural interpre- 
tation of the spin incommensurability in terms of antifer- 
romagnetic domains separated by charge stripes oriented 
along the principal axisi^2iiMi However, at lower doping, 
for 0.02 < x < 0.05, in the so-called spin-glass regime, 
where no superconductivity is observed, only two diag- 
onal peaks have been measured^ similar to the ones 
represented in Fig. 1. Even though these peaks could 
still arise from diagonal charge stripe formation, several 
arguments suggest that a spiral picture in this regime is 
more likely, as discussed in Refs. Illl36l38l 

These observations stimulated further investigations 
on the microscopic derivation of the 50(3) NLcr model 
(J2J), which can allow for the determination of the var- 
ious microscopic parameters. Klcc and Muramatsu 16 
considered the continuum field theory arising from a mi- 
croscopic spin-fermion model, where itinerant electrons 
are coupled via a Kondo-like interaction to localized 
spins, described by the Heisenberg model. The spin- 
fluctuations around the spiral configuration are in- 
cluded by allowing the vectors ni and ri2 to vary slowly 
on the lattice scale, and by adding a small ferromagnetic 
frustration L, which is afterward integrated out, 

S n + aL 

S ~ Vl + 2an • L + a 2 L 2 ' 

By also taking into account the coupling to the fcrmions 
and integrating them out, Klee and Muramatsu derived 
an effective action for the spin field which is the SO (3) 
NLer model © with an additional term 

S KM = J dtd 2 x[ Kk (d t n k ) 2 -p ka (8 a n k ) 2 } 

- J dtd 2 yis a (m • d a n 2 ) . (84) 

In the above action both the exchange J between the 
spins and the fermionic susceptibilities contribute to 
the coupling constants k and p a . If only the Heisen- 
berg interaction between the spins is considered, then 
«! = k 2 = K = l/{8Ja 2 [2 + J2 a cos(Q a a)}}, k 3 ~ 0, 
Plot = P2a = (J5 ,2 /4)cos(Q Q a), p 3a ~ and s a = 
(JS 2 /a) sm(Q a a). The subtle interplay between holes 
and spins is represented in the last term of the micro- 
scopically derived effective model (|84H . Indeed, since it 
is not positive definite, the weight of some field config- 
urations in the path integral will tend to infinity, hence 
leading to instabilities. In order to ensure the action to 
be at a minimum, one should impose the condition that 
the full coefficients s a — \a — {JS 2 /a) sin(Q Q a), where 
Xa is the holes' contribution, must vanish. This stabil- 
ity argument determines the spiral incommensurability 
Q a as a function of the microscopic parameters and the 
doping concentration^ 



A more general derivation of the stability condition has 
been proposed recently by Hasselmann et al. 11 In this ap- 
proach, Eq. (fH^|) is considered as the continuum limit of 
the Heisenberg model alone and the effect of doping is in- 
cluded within a minimum coupling of the order parame- 
ter to a gauge field B Q representing the dipolar character 
of the hole state, already emphasized by Shraiman and 
Siggia4iiii& It is then shown that the stability condition 

p ka d a n k ■ [B a ] D x n fe + s a ni • d a n 2 = 0, 

where [BJd denotes the ordered fraction of dipoles, re- 
lates the incommensurate vector Q to the hole densityipi 
in agreement with neutron scattering measurements in 
lanthanum cuprates^ Moreover, the dipolar frustration 
described within this minimal-coupling scheme renormal- 
izes the bare coefficients p ka of Eq. (|81J> leading the sys- 
tem toward a stable fixed point where p\ a = p 2a = p 3a = 
which corresponds to the model @ that we con- 
sidered. In this picture we can also determine the pa- 
rameter k = JS 2 /4c 2 , where c = ./cjjcX ~ 2^/2 J Sa and 
(c||/c^) 2 = cos(Qa). 11 As a consequence, we can now 
apply our previous results to the spin-glass phase of lan- 
thanum cuprates. 

B. Inverse mobility in cuprates 

Until now we evaluated the inverse mobility of the de- 
fect without specifying how this quantity can be accessed 
experimentally. As we explained above, in lightly doped 
cuprates the holes act at the same time as source and 
stabilizing mechanism of the dipolar frustration. When 
topological defects are present in the spiral spin texture, 
one could expect that the holes sit on top of the defects 
(single vortex or vortex-antivortex pair) to minimize the 
frustration. Thus, the measured in-plane inverse mobil- 
ity of the holes would be described by Eq. (|fl??)l with the 
damping coefficient given by Eq. I|7fi|l. However, this sce- 
nario should apply only for temperatures above 150 K, 
because below this temperature the experiments signal 
charge localization. 

In the case of cuprates, the magnon temperature E c is 
the antiferromagnetic coupling J « 1200 K measured at 
zero doping. Actually, a lower value is expected, if one 
takes into account the renormalization of the spin-spin 
interaction due to the disorder introduced by hole dop- 
ing and quantum effects. The resulting inverse mobility 
as a function of temperature for the case of a vortex- 
antivortex pair is reported in Fig. 6 for several values of 
E c and a. Observe that using a = 3.8A, as appropriate 
for cuprates, one obtains /ip 1 = 0.46 Vs/cm 2 . Inspection 
of Fig. 6 shows that already at E c = 1000 K the overall 
variation of /i^ 1 between 150 K and 300 K is of order of 
0.05 Vs/cm 2 , as observed experimentally — Moreover, in 
the case of the vortex-antivortex pair, an upper limit for 
the application of Eq. (|81|l is given by the temperature 
T v of the vortex-antivortex unbinding. By estimating 
T v ~ itJS 2 /2, in analogy with the XY model, we find 
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that T v is of order of 400 K for the values of E c used in 
Fig. 6. 



E c =1 106 K 

E r =1000 K 
E„=800 K 



0.07 
0.06 
0.05 
™E 0.04 



E c =1100 K, a=1 
E.=1 1 00 K, a=3 
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E r =1000 K, a=5 






FIG. 6: (Color online) Contribution 181H of the motion of a 
vortex-antivortex pair to the inverse mobility of the holes in 
cuprates. 



FIG. 7: (Color online) Contribution l|810 of the motion of 
a single-vortex defect to the inverse mobility of the holes in 
cuprates. 



VII. CONCLUSIONS 



The inverse mobility of a single defect is shown in Fig. 
7. The overall variation of the inverse mobility turns out 
to be quantitatively smaller in this case. As we explained 
in the previous section, the main difference between the 
temperature dependence of the inverse mobility obtained 
with the single vortex or with the vortex-antivortex pair 
is that in the former case yT 1 oc T 2 always, while in 
the latter case /i^ 1 evolves towards a linear behavior 
at a crossover temperature which depends on E c and 
a. As we discussed in Sec. II, the presence of the two 
kinds of defects depends on their energy, which scales as 
E[&i v ] ~ ln£ for the single vortex and £ , [ v E , 2u] ~ hid 
for the pair, where £ is the correlation length and d the 
distance between vortices, respectively. In the absence 
of disorder and at low temperatures, one would expect £ 
to be finite, but still large enough to prevent the forma- 
tion of free defects below the crossover temperature T v , 
where pairs start to unbind piiLSSi^i 2 ^ However, in Ref. 
ITU it has been argued that disorder leads to a strong 
reduction of the correlation length, and thus single de- 
fects start to proliferate already at temperatures lower 
than T v . Comparison with recent resistivity data seems 
to support this conclusion. Indeed, studies performed 
by Ando et alSL for compounds in the spin-glass regime 
indicate that the second derivative of the in-plane resis- 
tivity with respect to the temperature is positive up to 
~ 300!^, implying that p a t, ~ T v with rj > 1. Based 
on previous experiments^ we interpreted the resistivity 
data for cuprates in terms of the dissipative motion of 
vortex-antivortex pairs only, which gives rise to a linear 
resistivity. 12 However, in the light of these new data, one 
could speculate about the coexistence of single-vortices 
and vortex-antivortex pairs, which would consequently 
lead to a power-law behavior of the resistivity with a 
more complicated exponent, expected to be larger than 
one. 



We study here the properties of frustrated Heisen- 
berg spin systems in which a noncollinear spin state is 
formed at low temperatures. In the long- wavelength 
limit, the system is described by the 50(3) NLer model, 
and several differences arise with respect to the usual 
0(3) NLcr model adopted to describe collinear spin 
states. In particular, vortex-like excitations play a cru- 
cial role in determining the finite-temperature critical 
behavior >ii We concentrated on the contribution of these 
topological defects to transport properties. Our ap- 
proach extends to a non-Abelian field theory the well- 
known collective-coordinate method employed previously 
to study the dissipative mechanism in one- and two- 
dimensional systems^ We consider two kinds of topo- 
logical defects: a single vortex and a vortex-antivortex 
pair. We show that the interaction between the defect 
and the spin waves is described by a particle coupled to 
a bath of harmonic oscillators. The scattering of the de- 
fect by the magnons leads to its dissipative motion. We 
integrated out the bath and calculated the mobility of 
the defect. Quite generally, its temperature dependence 
is determined by the thermal activation of the magnons, 
which vanishes at zero temperature and follows, at higher 
temperatures, a power law whose exponent depends on 
the type of defect. In particular, we find that it is lin- 
ear for the vortex-antivortex pair and quadratic for the 
single vortex. We apply the model to describe transport 
in lightly doped lanthanum cuprates. Several theoretical 
and experimental studies suggest that in these systems 
a spiral state is formed at low temperaturesiii*i2*2Z^2i 
Our results for the mobility indicate indeed that a pos- 
sible mechanism for transport in these materials, for 
150K< T <400K, could be the dissipative motion of an 
electrical charge attached to a single vortex or a vortex- 
antivortex topological defect. 

Although we have applied the model to the particular 
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case of lightly doped cuprates, the approach presented 
here is quite general, and can be employed to investi- 
gate the role of topological defects in any frustrated spin 
system described by the SO(3) NLcr model. As far as 
the spin-glass phase of cuprate superconductors is con- 
cerned, the incommensurate peaks, the value of the resis- 
tivity, as well as its linear temperature dependence and 
anisotropy, might be explained within both the spiral and 
the stripe model>ii*i2i2L22 A theoretical prediction that 
would discriminate between these two scenarios, as well 
as its experimental realization, are still missing. This 
issue is currently under investigation. 43 
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APPENDIX A: ENERGY OF A SINGLE VORTEX 
AND OF A VORTEX- ANTIVORTEX PAIR 

In this appendix we will calculate the energy of the 
topological defects in the SO(3) NLcr model. Static so- 
lutions of the model obey the Laplace equation 



V 2 *(r) = 0. 



(Al) 



A single- vortex solution centered at R = (A, Y) has the 
form 



*i t ,(r, R) = arctan 



x — X 

y-Y 



(A2) 



Its energy is given by 

E[9iv] = Mc 2 [ d 2 r(Wi„) 2 = 2^c 2 ln-. 

J a 

A bound vortex-antivortex pair described by 

*2v = *i„(r, Ri) - *i„(r, Ra), (A3) 

with 'I'i,, given by Eq. I|A2(1 . is also a solution of Eq. 
(|A1I) . The vortex-antivortex defect can be written in a 
more compact form as 



^2v = arctan 



[dx (r-R)] z 
(r-R) 2 -d 2 /4 



where the center of mass and relative coordinate, respec- 
tively, are given by 

1 



In order to evaluate the energy of the vortex-antivortex 
pair 

E[^ 2v ] = Nc 2 J d 2 r(W 2 „) 2 , 

we use Eq. I|A3|) . which yields 

E[* 2v ] = M<?{h x + J 22 - 2/ 12 ), (A4) 

where 

Ijji = J d 2 r(W (j) )(V* (/) ). 

and ^(i' 2 ) = »I< lt ,(r, R(! 2 )). It is easy to show that In 
and I22 are equal, 



In = I22 = 27rln -. 

a 



(A5) 



The integral 7i 2 is highly non-trivial. After some calcu- 
lations, it can be expressed in the form 



1 



1 1 



/|J == 2.1 d r (r-R!) 2 ' 2 J - >-R 2 ) 2 



r / d z v- 



1 



- i( Rl -R 2 ) 2 y d 2 r- 



1 



r-R!) 2 (r-R 2 



,(A6) 



The first two integrals on the right hand side of Eq. (|A6|) 

are identical and equal to In, whereas the last one must 

(3) 

be evaluated separately. Let us denote it as I{ 2 . After 
regularization, it becomes 



l[f = I d 2 r 



1 



[(r-R 1 ) 2 + a 2 ][(r-R 2 ) 2 + a 2 ]- 



By introducing new coordinates with r — R 2 — ► r, the 
last integral acquires the form 



ig> = d 2 r 



1 



[r 2 + a 2 l[(r + d) 2 + a 2 



(A7) 



In order to simplify Eq. I|A7|I . we introduce polar coordi- 
nates x — r cos ip,y — r sin ip 



r(3) 



dr- 



d^p 



+ a 2 J r 2 + 2rdcos(ip — </>) + d 2 + < 



where (d, <p) are polar coordinates of d. Integrating over 
the angle ip, we obtain 



2?r 



r 2 + a 2 ^J[r 2 - d 2 ) 2 + 2a 2 {r 2 + d 2 ) + a 4 ' 



After substituting u = l/(r 2 + a 2 ), the above integral 
reads 



r(3) _ _ 



l/(f 2 +a 2 ) 



d:U 
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Straightforward calculations give 



which allows us to write 



r(3) _ 
J 12 — 



In (A k + V^" 

dVd 2 + 2a 2 I 4a 4 L 



d + v 7 ^ 2 + 2a 2 } 



4tt 



dVd 2 + 2a 2 a 
Substitution of Eq. l|A8j) into Eq. I|A6|I yields 



ln-.(A8) 



hi = 2tt 



In - In — 

a Vd 2 + 2a 2 a 



(A9) 



After inserting the last relation together with Eq. (|A5|) 
into Eq. I|A4(1 . we finally obtain the energy of the vortex- 
antivortex pair 

25[* a „] = — ^L=AAc 2 In - « 4^c 2 In - , (A10) 
VgP + 2a 2 a a 

which shows that the energy of the defect pair is finite. 



A« = -tr^ - a ' — »™ " » - b 



>1 - 4"M Cr I COS— -ZTO-Crsm— | rr 



2 



2 



cos + im ■ trsin ) \ m b d^ v ( 1 



— tr < cr I cos — ZTO • cr sin | n 



cos — - + im • asm — - I (a ■ e)> m b d^ v 



-tv ^ g [a ■ e) [ cos — im ■ a sm — 



d b ( cos + im ■ ersin ) ^ m b d^ v 



1 . / „~ - / *t> 

— tr < a (a ■ e) cos im • <r sm — 

16 \ v ' V 2 2 



a ( cos — — h im ■ asm ) (a ■ e)f m d^ v . 



APPENDIX B: DYNAMICS OF THE 
FLUCTUATIONS AROUND THE DEFECT 

Using the identities 



After some algebra, the last equation acquires the form 
- \m a d^ v (l - j) + \e abc e h m c d^ v 



a a 



5" 



b , ■ abc c 

ie a , 



/ I \ a a ■ a a 

exp —a - a) — cos — + i sin — . 

\2 J 2 a 2 

where a = \a\, together with Eq. one can show that 



d^9s = -m a a a g s d^ v , 



(Bl) 



i a ■ e 



2 e 



id" 

E 



Here e stands for |e|. By inserting Eqs. IjBlfl into the field 
A® given by 10), we obtain 



A a = A a 4- A a 4- A a 4- A a 



(B2) 



with 



A a 



A a 



A a 



A a 



1 

4 
1 

2s 



ti{d a g- 1 g- 1 d b g s g e )m%* v , 



tr(d a g^d^e" sin-, 



-ti(d a g- 1 d b ) y ^ ' sin 



2' 
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tv(d a g- i d b d hl eg e )e b 



e a e b m b d^ v - -e 2 m a d^ v + 0(e J ). (B3) 



Analogously, one can show that 
1 



, ^ 0(E* 



A a — 



e a (e ■ d^e) 



(B4) 



By substituting Eqs. l|B3j) and (|B4|) into the field A a 
given by l|B2|) , and retaining terms up to the second order 
in e, we obtain 

Al = l -m a d^ v f 1 - y) + \e a e b m b d^ v 

+ l -e abc e b m c d^ v + ^d,s a + ie afcc e h ^ £ c ,(B5) 



and 



A*Al = l( d ^ v f + 1(8^ + ±m a d^ a d^ v 



(B6) 



The parameter e is small, e <g; 1, because g e describes 
fluctuations around the defect. Using the properties of 
the Pauli matrices, we find 

& = 1 + |<?- £-i £ 2 + 0( £ 3 ), 



APPENDIX C: EVALUATION OF THE KERNEL 

In this appendix we will express the kernel K defined 
by Eq. I|4l)|l as a functional integral. First, we divide the 
time interval [0, t] into (m — 1) subintervals of length e, 
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so t = [m — l)e, and use (to — 1) completeness relations 
between the (to — 1) exponential functions, 

K{x,a*;y,0;t) = (xl, a\e~ ^ \y , 0) 



= / d Xm-i... / d 2 x 1 



d 2 Cm-i f d 2 Ci 



T 2N 



T 2N 



j 7r^ 

(XmCml^- h |^m— lCm— 1/ 

(x TO _i(f m _i|e _ ~|x m _2Cm-2)---(xiCi|e"~ s ~ |xoCo), 



where 



x m = x, x = y 

Cm = a, Co = 



(CI) 



By inserting m completeness relations in the momen- 
tum representation into Eq. (|C1() . we obtain 



f m — 1 



m 

[J{x fe Cfe|e-^|P fe Cfe-i)(Pfc|x fc _i). (C2) 



fe=i 



The matrix element (x/cCfe | e n |PfcCfe— 1)> (A; = 1, m) , 
can be evaluated using that e<t. It reads 



(xfcCfc|e '^IPfeCfe-i) = (xfc|Pfe>(Cfc|Cfe-i) 



exp 



-i/(x fc ,P fc ;C fc *,Cfc-i) 



(C3) 



where 



Ut T> t* t \— ( X feCfc|-ff|PfcCfc-l) 



(xfclPfcXaiC 



k-1) 



Substituting the matrix element given by Eq. (|C3|) into 
Eq. JU2J), recalling that 



(Plx) 



1 I 

exp I — — P • x 



2nh 



and using the properties of the overlap of coherent states, 
we find that the kernel acquires the form 



K(x,d?;y,/?;t) 



d 2 V 

2-Kh 



n 



ri 2 X ; rf 2 P, ri 2 Q 

2nh tt 2N 



{m ^ 
fe=i 

^ [P fc (x fc - xfc_i) - efffxfc, P fe ; C, Cfc-i) 



ent state representation for the bath of harmonic oscilla- 
tors, we find 

-» - P 2 

H(x fe ,P fe ;Cfe,Cfc-l) = ^TT + E ^""Cnm.iC^i 



mn,'i—a,b 



Pfe 

M 



ETd z^- 1 _ D Afe-i Afe* 



which, after insertion into Eq. (|C4|) and integration over 
the momenta Pfe, yields 



K(x,a*;y,M 



M 



ill / m— 1 ,2 j2 / 



{m 
fe=i 



with 



pq.b 



Now, we consider the continuum limit, e — > 0, of the last 
equation. Using the boundary conditions (|C1|1 . the first 
term in the exponential in (|C5|) becomes 



ex p {e \ • (3 - Cfe-i) - a • (d - d-i)] I 



fc=0 



fc=l 



= exp \ -[Co • (cr - Co) - c* • (c™ - Cm-i)] 



cx p { ^ E tc* • (Ck+i -a)- a- (d - | 

exp ^i|/3| 2 - i|a| 2 ) exp • C*(0) + \d* ■ C(i) 



exp<ji / d*'(C • C* - C* • C) 



(C4) 



In order to integrate over the momenta Pfe, we have to 
explicitly calculate H(x.k, pfej Q, Cfe-i)- Using the coher- 



The other terms in Eq. (|C5|I can be trivially written in 
the continuum limit, yielding then Eq. 1)41 [I. 



APPENDIX D: INITIAL DENSITY MATRIX FOR 
THE BATH IN THE COHERENT STATE 
REPRESENTATION 

It remains to evaluate the matrix elements of the initial 
density matrix for the bath in the coherent state repre- 
sentation 



p B 0*J',O) = 0\p B (O)0') 
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with 



APPENDIX E: EQUATIONS OF MOTION 



PB{0) = -\{e 



-Uuj pq (al „a p „+£>+ b pq ) 



where the partition function reads 



Z = ti 



~UT, pq ^pq(a-l q a. Pq +bl q b pq ) 



n i 



V! 



Since the baths a and b are not coupled, the total density 
matrix is the product 



PB 



0*J , ,O)= P B,a0:Ja,O)PB,b0* b J'b,O) 



of the density matrices pB,a and pB,a for the baths a and 
b, respectively, given by 



l"l 



for the bath a, and analogously for the bath b, with the 
partition function 



Zb = Z a = tr 



""Ep, u pq {al q a pq )\ =TT ( 1 



We shall evaluate the previous matrix element by insert- 
ing two unity operators in the occupation number rep- 
resentation for the bath a (in the following, to simplify 
notation we omit the index a in the /3, f3' and n) 

PB , a 0*J',O) = ±l[(f3 pg \e- u ^ & U & ^' pq ) (Dl) 



p<-i 



rU E <* 



Z a 



'pq | Tlpq 



{n m \e- u ^ & U^\ n ' pq ){n' \[3' pq >. 



Now, we use the scalar product of the states which de- 
fine the coherent state and occupation number represen- 
tations 



(A 



pq\' l pql 



-e : 



iipq . 



Substituting the above expression into Eq. IjDlfl . we ob- 
tain 



nra \ 

y— 



P9 \ 



I | I /^pg I 

Z a J - L \ ~~ 2 2~ 

exp [/3* 9 /3p g exp(-L/w pg )] . 



(D2) 



Our next step is to solve equations of motion l|47|) . In 
order to achieve this aim, we introduce the ansatz 

Cnm(r) = e~ iu '™ r \P nm + W nmM^)hl], 

kl^nm 

C* ( T ) = e luJnmT \a* e ~ iUnmt 

+ W^uWe-^'afa], (El) 

kl^nrn 

where the functionals W and W will be determined from 
the equations of motion and < r < t. By substituting 
the first time derivative of Eqs. (|E1|) into the equations of 
motion (|47() . we obtain the expressions which determine 
the time evolution of the functionals W and W 



W nm ,kl - W° m kl + ^ W nm,pq W pq,kh 
pq 



W nm ,kl = ~W° m<kl - ^ Wn m ,p q Wp qt k 



with 



WL,«fcr] = iiD^.we-*^""-""^. (E2) 

Notice that = -W°^ H . Because W and W 

must satisfy the boundary conditions W(x, 0) = and 
W(x, t) — 0, we have 

W nm , kl [x,r] = [ dt'Wl mM [*A 
Jo 

+ E / dt'W^ m<pq [^t']W pq , kl [^t'}, 



W'nm.wlx.r] = / dt'W°[x,t'} 



+ E / dt 'wL,p q [x,t']w pq , kl [*,t'], 

pq Jt 

(E3) 

which, using the Born approximation, acquires the form 

Wnm,«[x,r] = f dt'W% m>kl [x,t'} 

Jo 

+ E / MwLmM i dt"W° qM {x,t'}, 
,„„ Jo Jo 



W nm , k i[^r} = j dt'W% mM [x,t'} 
+ E / ^Vps [x, f ] jf *"W&,w [x, ■ (E4) 
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The functions 7 appearing in Eq. (|44|l obey the equations with 
of motion l|47|l with the boundary conditions (|46[) . We 
solve them by introducing the ansatz 



^ nm.kl ^^nm.kh vv nm,kl yv nm,kl' 



-iu nm r r iuj nrn t 



W° ur = W°* 
y r nm,kl r * n 

The boundary values of the functionals obey the relations 
WW;[x,0] = W* H [x,0], 



7«m(T) = 

+ E W nmM {T)a kl e l ^% 

kl^nm 

with the conditions = and W(0) = 0. By in- which will be used later. From Eqs. (|E1(I and 1|E5|) the 

serting this ansatz into the corresponding equations of boundary terms read 



(E7) 



motion, we find, after some algebra, 
W nmM [x,T] = / dtW nmM [x,t] 



E / dt'W^ pq [x,t']W pq ,kl 



(nm(t) = Primed™* + E W nm,kl (x, t^"™* kl , 

C m (0) = < m e~ iaW + E^»™,^(x,0)e— 
7nm(0) = ewe*"-"* +E^™^«(y,0)e ia " !!t afci, 
7™(*) = t^"' + E ^n m ,«(y, tje**"-'^. (E8) 



(E6) 



J 



APPENDIX F: EVALUATION OF T 
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In this appendix we evaluate the diagonal elements of the matrix T = T a + T b , where the elements T^ m nm are 
given by 



r . = ir 

nm,nm ^ 



W„ m ,„ m [x,t] + W, 



nm.nm [■**-) 



0] + ^„*™,„ ro [y 5 0] + W* mtnm [y, t}] + 1 E [w^k 0] + W M ,nm[x, i] 



W^* [v 01 + W* \v t] 



and the ones of matrix r b are obtained from the latter by the substitution D, 



.kl 



D 



kl.nm — ^nmM 



(Fl) 



Using the 



Born approximation for the functionals W and W given by Eq. I]E4[I. and the form of the functionals W and VF°, 
defined after Eq. IjEll) . we find 



nm,pq 



j(u nm -OJ pq )(t'~t") 



W nm , nm [x,t] = - E / W dt"9(t' -t")x^(t')^(t")D^ lpq D- m 
VF„ m ,„ m [x,0] = - E I dt' ( dt"6{t" -t')x»{t')x v {t")D» m>p ^ 



p,,v,pq 



(F2) 



Using Eq. I|E3|) . as well as its complex conjugate evaluated at y and retaining only the terms quadratic in the coupling 
constants, we can write the last term in Eq. I|F1|) as 



\ E^™(x, 0) + W„(x, t)][W: m , pq (y, 0) + W£, nm (y, t)] ^-j dt 1 [ dt" \w° mtPq [K, t'] + , nm [x, t'] 



JO 



L J .. Jo Jo 



-uw)(t'-t") 
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= J df dt"6(t' - t")Dl m vq D v n * m vq ^( t ')^( f ") e i(^w-^ m )(t'-i") + ^(f^^ e -i(w pq -u nm )(t'-t'')^ (F3) 

Substituting Eqs. (|F2|) and (|F3|) into Eq. |IF1|) . we obtain the diagonal elements of the matrix T a , in the lowest 
order in the coupling constants 



y V / nm,pq nm,pq^ 

yielding the diagonal elements of the matrix T 



r nm , nm = - £ / dt' [ ^(t'-Ofe.M^^ + ^M^^)^^)-^^)] 

,, „ ™ Jo Jo 



[i"(t") 



,»Km-Ujj)(t'-*") _ y v (t")p~ K u nm.-U pg )(t'-t") 



(F4) 



(F5) 



APPENDIX G: EVALUATION OF THE COUPLING CONSTANTS 

In this appendix we calculate the coupling constants 

where the wave functions are given by 



in I'd 



(m > 0). 



By expressing the gradient operator in polar coordinates we find that 

p oc 

Gfcmfc'i = 7r« TO _i,i / drr[i? i (fcV)5 m (fcr)-mi 1 , i (fcV) J F 1 m (fcr)] 
Jo 

poo 

+ nS m -i,-i / drr[Ri(k'r)S m (kr) + mRi(k'r)F m (kr)], 
Jo 

/>oo 

Ct'i - «r<y m -J,l / drr[i? i (fc'r)5 m (fcr)+TOi? ; (fcV)F m (fcr)] 



wr£m-i,-i / drr[Ri(k'r)S m (kr) - mRi(k'r)F m (kr)], 
Jo 



where 



Sj (At) 



fe 
2/ 

* /A 
2V 2l 

,-2i6, rr(2) 



^[^(ArJ + e 2 ^ 15 ^)], 



mr) = ^vS [ ^ (1)(fcr)+e_2Wijff ' (2)(fcr)] - 



In evaluating the expression (|G3(1 one can use the asymptotic form of the Hankel functions 

Hg>(kr) = 
H%\kr)- 



nkr 



g ifer e -if (m+l/2) 



nkr 



e -ikr £ i^ (m+1/2)^ 



(Gl) 



(G2) 



(G3) 

(G4) 

(G5) 
(G6) 
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However, the terms coming from F m (kr) would then be divergent at r = 0. This divergence is an artifact of 
approximating up to r — the true solution of the scattering problem i|2U[l with the functions (|G2fl . At small r indeed 
a better approximation for the radial part of the functions rjkm is provided by the Bessel functions J m (kr), which are 
regular at r = 0. Then wc find that 



drrH^ 2 \k'r)H^ 2 \kr) = 
drrHjP(k'r)H^(kr) = 



2 e ±i*( t i+v+l)/2 



ir\/kk' 
nvkk' 



irS(k + k')±V 



k + k' 



nS(k -k')+V 



1 



k-k' 



um i 

drJ m+1 (k'r)J m (kr) = JJJy^I^k' - k) + -S(k - k'). 



(G7) 



Because the terms in S(u>, u>') proportional to S(lu±lo') do not contribute to the damping matrix Ijtjlfl . we can directly 
discard them from the definition of the coupling constants, which are then given by 



k 



G km,k'l ~ J$l,m±l | i^hlV k + k , - A hn V k _ jj 



TO7TV kk' 



(k') m 1 



mirykk' 



^ ri e(k'-k) + -s(k>-k) 



L (l) r 1 

lm k + k' 



2k 

A^y 1 
lm k-k' 



At 



(k>) 



miry'kk' 



^n-e(fc - k') + —5(k - k') 



e { k>-k) + -8{k>-k) 



(G8) 



where 



a(1) _ An{l+m)/2 p -2id m , (> -iir(l+m)/2 e> 2i5i a (2) _ ix(l-m)/2 , -iir(l-m)/2 2i(8i-S m ) 
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